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ABSTRACT 

We show that the exoplanet HAT-P-7b has an extremely tilted orbit, with a true angle of at least 86° with 
respect to its parent star's equatorial plane, and a strong possibility of retrograde motion. We also report 
evidence for an additional planet or companion star. The evidence for the unparalleled orbit and the third body 
is based on precise observations of the star's apparent radial velocity. The anomalous radial velocity due to 
rotation (the Rossiter-McLaughlin effect) was found to be a blueshift during the first half of the transit and 
a redshift during the second half, an inversion of the usual pattern, implying that the angle between the sky- 
projected orbital and stellar angular momentum vectors is 1 82 .°5 ± 9 .°4. The third body is implicated by excess 
radial-velocity variation of the host star over 2 yr. Some possible explanations for the tilted orbit are a close 
encounter with another planet, the Kozai effect, and resonant capture by an inward-migrating outer planet. 

Subject headings: planetary systems — planetary systems: formation — stars: individual (HAT-P-7) — 
stars: rotation 



1. INTRODUCTION 

In the Solar system, the planetary orbits are well-aligned 
and prograde, revolving in the same direction as the rotation 
of the Sun. This fact inspired the "nebular hypothesis" that the 
Sun and planets formed from a single spinning disk (Laplace 
1796). One might also expect exoplanetary orbits to be well- 
aligned with their parent stars, and indeed this is true of most 
systems for which it has been possible to compare the direc- 
tions of orbital motion and stellar rotation (Fabrycky & Winn 
2009, Le Bouquin et al. 2009). However, there are at least 3 
exoplanets for which the orbit is tilted by a larger angle than 
any of the planets in the Solar system: XO-3b (Hebrard et 
al. 2008, Winn et al. 2009a), HD 80606b (Moutou et al. 2009, 
Pont et al. 2009, Winn et al. 2009b), and WASP- 14b (Johnson 
et al. 2009). 

Still, all of those systems are consistent with prograde or- 
bits, with the largest minimum angle between the stellar- 
rotational and orbital angular momentum vectors of about 
37°, for XO-3b (Winn et al. 2009a). The reason why only the 
minimum angle is known is that the evidence for misalign- 
ment is based on the eponymous effect of Rossiter (1924) and 
McLaughlin (1924), an anomalous Doppler shift observed 
during planetary transits that is sensitive only to the angle be- 
tween the sky projections of the two vectors. The true spin- 
orbit angle may be larger, depending on the unknown inclina- 
tion angle of the stellar rotation axis with respect to the line 
of sight. 

In this Letter we present evidence of a very large spin-orbit 
misalignment for HAT-P-7b, a planet of mass 1.8 Mj up and 
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radius 1.4 Rj llp in a 2.2-day orbit around an F6V star with 
mass 1.5 M Q and radius 1.8 R (Pal et al. 2008). We find the 
angle between the sky-projected angular momentum vectors 
to be 182.°5 ± 9.°4. Furthermore we show that the true angle 
ip between those vectors is likely greater than 86°, indicating 
that the orbit is either retrograde (ip > 90°) or nearly polar 
(tp w 90°). We also present evidence for a third body in the 
system, which may be an additional planet or a companion 
star. We present spectroscopic data in § 2, photometric data 
in § 3, a joint analysis of both types of data in § 4, and a 
discussion of the results in § 5. 

2. RADIAL VELOCITIES 

We observed HAT-P-7 with the High Resolution Spectro- 
graph (HIRES) on the Keck I 10m telescope, and the High 
Dispersion Spectrograph (HDS) on the Subaru 8m telescope. 
The planet's discoverers (Pal et al. 2008; hereafter, P08) ob- 
tained 8 HIRES spectra in 2007, to which we add 9 spectra 
from 2009. All but one of the HIRES spectra were acquired 
outside of transits. Of the 49 HDS spectra, 9 were obtained 
on 2009 June 17 and 40 were obtained on 2009 July 1. The 
second of these nights spanned a transit. 

The instrument settings and observing procedures in both 
2007 and 2009 were identical to those used by the California 
Planet Search (CPS; Howard et al. 2009). We placed an io- 
dine gas absorption cell into the optical path, to calibrate the 
instrumental response and wavelength scale. The radial ve- 
locity (RV) of each spectrum was measured with respect to an 
iodine-free template spectrum, using the algorithm of Butler 
et al. (2006) with subsequent improvements. Measurement 
errors were estimated from the scatter in the fits to individual 
spectral segments spanning a few Angstroms. The RVs are 
given in Table 1 . 

2.1. Evidence for a third body 

Fig. 1 shows the RVs over the 2 yr span of the observations. 
Fig. 2 shows the RVs as a function of orbital phase, fitted 
with 2 different models. The first model is a single Keplerian 
orbit, representing the signal of the known planet. The sec- 
ond model has an additional parameter 7 representing an ex- 
tra radial acceleration. The second model gives a better fit to 
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Fig. l.— Long-term radial velocity variation of HAT-P-7. (a) Measured RVs. (b) Residuals (observed - calculated) be- 
tween the data and the best-fitting single-planet model. Light blue and dark blue points are HIRES data from 2007 and 2009, 
respectively. 



the data, with a root-mean-squared (rms) residual of 7 m s -1 
as compared to 21 m s" 1 for the first model. The RVs from 
2009 are systematically redshifted by approximately 40 ms" 1 
compared to RVs from 2007, as evident from the residuals 
shown in Figs. 1(b) and 2(b). This shift is highly significant, 
as the CPS has demonstrated a long-term stability of 2 m s -1 
or better using HIRES and the same reduction codes used here 
(Howard et al. 2009). 

This RV trend is evidence for an additional companion. 
Given the limited time coverage of our observations (two clus- 
ters of points separated by 2 yr), the data are compatible with 
nearly any period longer than a few months. A constant ac- 
celeration is the simplest model that fits the excess RV vari- 
ability, and under that assumption we may give an order-of- 
magnitude relation relating 7 to some properties of the com- 
panion, 

M f sini c ^ 2 = ((U2 l ±0.014)M Jup AIT 2 , (1) 

a* G 

where M c is the companion mass, i c its orbital inclination rel- 
ative to the line of sight, a c its orbital distance, and the numer- 
ical value is based on our model-fitting results (see § 4). 

2.2. Evidence for a spin-orbit misalignment 

Fig. 3(a) shows the RV data spanning the transit, after sub- 
tracting the orbital RV as computed with the best-fitting model 
including 7. We interpret the "anomalous" RV variation dur- 
ing the transit as the Rossiter-McLaughlin (RM) effect, the 
asymmetry in the spectral lines due to the partial eclipse of the 
rotating photosphere. In the context of eclipsing binary stars, 
the RM effect was predicted by Holt (1893) and observed 
definitively by Rossiter (1924) and McLaughlin (1924). For 
exoplanets, the RM effect was first observed by Queloz et 
al. (2000), and its use in assessing spin-orbit alignment has 
been expounded by Ohta et al. (2005) and Gaudi & Winn 
(2007). 

A transiting planet in a well-aligned prograde orbit would 
first pass in front of the blueshifted (approaching) half of the 



star, causing an anomalous redshift of the observed starlight. 
Then, the planet would cross to the redshifted (receding) 
half of the star, causing an anomalous blueshift. In contrast, 
Fig. 3(a) shows a blueshift followed by a redshift: an inver- 
sion of the effect just described. We may conclude, even with- 
out any modeling, that the orbital "north pole" and the stellar 
"north pole" point in nearly opposite directions on the sky. 

3. PHOTOMETRY 

For a quantitative analysis of the RM effect we wanted to 
model both the photometric and spectroscopic transit signals. 
For this purpose we supplemented the RV data with the most 
precise transit light curve available to us, shown in Fig. 3(c). 
This light curve is based on observations on UT 2008 Sep 22 
in the Sloan i bandpass, with the Fred L. Whipple 1.2m tele- 
scope and Keplercam detector, under the auspices of the Tran- 
sit Light Curve project (Holman et al. 2006, Winn et al. 2007). 

Reduction of the CCD images involved standard procedures 
for bias subtraction and flat-field division. Differential aper- 
ture photometry was performed for HAT-P-7 and 7 compari- 
son stars. No evidence was found for time-correlated noise us- 
ing the "time-averaging" method of Pont et al. (2006), as im- 
plemented by Winn et al. (2009c). The data shown in Fig. 3(c) 
were corrected for differential extinction as explained in § 4. 

4. JOINT ANALYSIS 

We fitted a model to the photometric and RV data in or- 
der to derive quantitative constraints on the angle A between 
the sky projections of the orbital and stellar-rotational angu- 
lar momentum vectors. This angle is defined such that A = 0° 
when the sky -projected vectors are parallel and A = 180° when 
they are antiparallel. Our model for the RM effect was based 
on the technique of Winn et al. (2005): we simulated spectra 
exhibiting the RM effect at various transit phases, and then 
measured the apparent RV of the simulated spectra using the 
same Doppler code that is used on actual data. This allowed 
us to relate the anomalous RV to the parameters and positions 
of the star and planet. 
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Fig. 2. — Phased radial velocity variation of HAT-P-7. (a) Assuming a single Keplerian orbit, (b) Residuals, (c) With an extra 
parameter 7 representing a constant radial acceleration, (d) Residuals. The circles are HIRES data (light blue from 2007, dark 
blue from 2009), the green triangles are HDS data from 2009 June 17, and the red squares are HDS data from 2009 July 01. 
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The RV model was the sum of the Keplerian RV and the 
anomalous RV due to the RM effect. The photometric model 
was based on the analytic equation for the flux of a quadrati- 
cally limb-darkened disk with a circular obstruction (Mandel 
& Agol 2002). As a compromise between fixing the limb- 
darkening coefficients U\ and 112 at theoretically calculated 
values, and giving them complete freedom, we fixed U\—Ui 
at the tabulated value of 0.3846 (Claret 2004) and allowed 
u\+U2 to be a free parameter. We also included a free pa- 
rameter for the coefficient of differential airmass extinction 
between HAT-P-7 and the ensemble of comparison stars. 

We determined the best values of the model parameters and 
their 68.3% confidence limits using a Markov Chain Monte 
Carlo algorithm, as described in our previous works (see, e.g., 
Winn et al. 2009a). The likelihood function was given by 



N f r n 2 

/i(obs)-/;(calc) 



E 



N r -|2 

V|(obs) — v/(calc) 



E 



cr v j 



(2) 

in a self-explanatory notation, with er/., chosen to be 0.00136, 
and <7 V ,; chosen to be the quadrature sum of the RV measure- 
ment error and a "stellar jitter" term of 9.3 m s . These 
choices led to x 2 = Ndof for the minimum-x 2 model. A Gaus- 
sian prior constraint was imposed upon the orbital period 
based on the precise measurement of P08. 

Table 2 gives the results for the model parameters. In par- 
ticular, the result for A is 182.5 ±9.4 deg, close to antiparallel, 
as anticipated from the qualitative discussion of § 2. 



5. DISCUSSION 
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fig. 3.— The spectroscopic and photometric transit of HAT-P-7b. (a) The anomalous RV, defined as the output of the 
Doppler code minus the orbital RV. We observed a blueshift in the first half of the transit, and a redshift in the second half of the 
transit, demonstrating that the sky projections of the orbital and stellar angular momentum vectors point in opposite directions, 
(b) Residuals. Red squares are HDS data from 2009 July 1, and blue circles are HIRES data obtained on various nights in 2007 
and 2009. (c) The relative flux, observed in the Sloan i band with the FLWO 1.2m telescope and Keplercam. (d) Residuals. In 
panels (a) and (b), the gray line shows the best fitting model. 



Our finding for A is strongly suggestive of retrograde mo- 
tion, in which the orbital motion and stellar rotation are in 
opposite directions. However, it must be remembered that A 
refers to the angle between the sky-projected angular momen- 
tum vectors. The true angle ip between the vectors is given by 



cosip = cos /* cos /+ sin/* sin /cos A, 



(3) 



where i and /* are the line-of-sight inclinations of the or- 
bital and stellar angular momentum vectors, respectively. Al- 
though i is known precisely from the transit data, /* is un- 
known. 

Supposing /* to be drawn from an "isotropic" distribu- 
tion (uniform in cos/*), the data demand that ip > 86.°3 with 
99.73% confidence. Thus, under this assumption, a retrograde 



orbit is strongly favored, although a nearly-polar and barely- 
prograde orbit cannot be ruled out. 

In fact there is circumstantial evidence that /* is small and 
consequently the orbit of HAT-P-7b is nearly polar (ip « 
90°). The star's projected rotation rate is unusually low for 
such a hot star: vsin/* = 4.9^9 km s" 1 in our model, or 
3.8 ±0.5 km s" 1 based on the line profile analysis of P08; 
and r e ff = 6350 ± 80 K according to P08. In the SPOCS 
catalog of dwarf stars with well-determined spectroscopic 
properties (Valenti & Fischer 2005), only 2 of 37 stars with 
T eS = 6350 ± 100 K have vsin/* < 4.9 km s" 1 . 

Based on this catalog, the mean rotation rate v for such hot 
stars is about 15 km s" 1 . As an alternate approach to con- 
straining ip, we assumed the rotation velocity v is drawn from 
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a Gaussian distribution with mean 15 km s -1 and standard de- 
viation 3 km s" 1 . The result is i\> = 94. 6^3 q deg with 68.3% 
confidence, and i/j > 86.° 1 with 99.73% confidence. This anal- 
ysis favors nearly-polar and retrograde orbits. However, one 
wonders whether HAT-P-7 should be expected to have a "typi- 
cal" rotation rate, given the existence of its short-period planet 
on a bizarre orbit. Another caveat is that we found the scaled 
semimajor axis a/R+ to be about ltr smaller than the find- 
ing of P08, suggesting the star is somewhat larger and more 
evolved, which would correspond to a slower expected rota- 
tion rate. 

Determining i+ directly may be possible by measuring 
and interpreting asteroseismological oscillations (Gizon & 
Solanki 2003), or photometric modulations produced by 
starspots (see, e.g., Henry & Winn 2008). By good fortune, 
HAT-P-7 is in the field of view of the Kepler satellite, which 
is capable of precise long-term photometry and may be able 
to accomplish these tasks (Borucki et al. 2009). 

The extraordinary orbit of HAT-P-7b presents an extreme 
case for theories of planet formation and subsequent orbital 
evolution. HAT-P-7b is a "hot Jupiter" and presumably mi- 
grated inward toward the star after its formation. A prevail- 
ing migration theory involves tidal interactions with the proto- 
planetary disk, but such interactions would probably not per- 
turb the initial coplanarity of the system, and might even bring 
the system into closer alignment (Lubow & Ogilvie 2001, 
Cresswell et al. 2007). More promising to explain HAT-P-7b 
are scenarios involving few-body dynamics, as those scenar- 
ios are expected to produce misalignments. In one scenario, 
close encounters between planets throw a planet inward, 
where its orbit is ultimately shrunk and circularized by tidal 
dissipation (Chatterjee et al. 2008, Juric & Tremaine 2008). 
Another idea is based on the Kozai (1962) effect, whereby the 
gravitational force from a distant body on a highly inclined 
orbit strongly modulates an inner planet's orbital eccentricity 
and inclination (Fabrycky & Tremaine 2007). Recent calcu- 
lations showed that a combination of planet-planet scattering, 
the Kozai effect, and tidal friction can lead to nearly-circular 
retrograde orbits (Nagasawa et al. 2008). A third proposed 
scenario involves an inward-migrating outer planet that cap- 
tures an inner planet into a mean motion resonance; if the 
inner planet avoids being ejected or consumed by the star, it 
may be released on a nearly-circular retrograde orbit (Yu & 
Tremaine 2001). 

The prospect of explaining HAT-P-7b's orbit through few- 
body dynamics lends extra importance to measuring the mass 



and orbital parameters of the third body. If it turns out to 
be a planet, then HAT-P-7b will be only the second known 
case of a transiting planet accompanied by another planet, 
the first being HAT-P-13b (Bakos et al. 2009). Such systems 
are highly desirable because the unusually precise measure- 
ments enabled by transit observations can be used to deter- 
mine whether the orbits are coplanar and give clues about the 
system's dynamical history (Fabrycky 2009). 

Note added after submission. — Narita et al. (2009) report 
independent evidence for a retrograde or polar orbit of HAT- 
P-7b, based on Subaru/HDS spectra spanning the transit of 
2008 May 30. 
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TABLE 1 

Relative Radial Velocity Measurements of HAT-P-7 



HJD 


RV[m s~'] 


Error [m s 1 ] 


Spec. 11 


2454336.73960 


111.08 


1.72 




2454336.85367 


58.89 


1.78 




2454337.76212 


-236.06 


1.70 


1 


2454338.77440 


151.06 


1.54 




2454338.85456 


131.12 


1.57 




2454339.89886 


-256.42 


1.83 




2454343.83180 


-162.19 


1.75 




2454344.98805 


84.42 


2.15 




2454983.99020 


-77.44 


2.20 


1 


2454985.02095 


110.41 


1.96 


1 


2454987.01053 


189.80 


2.23 


1 


2454988.00940 


-221.82 


2.28 




2455014.95322 


-41.51 


2.08 


1 


2455016.05508 


4.10 


2.23 




2455042.03696 


189.91 


2.53 




2455042.98594 


-228.74 


2.40 


' 


2455044.07369 


191.34 


2.55 


1 


2455000.03550 


196.69 


12.26 


2 


2455000.07080 


184.21 


11.29 


2 


2455000.07507 


210.15 


11.25 


2 


2455000.07933 


195.95 


11.94 


2 


2455000.08361 


193.11 


13.30 


2 


2455000.08787 


225.66 


12.45 


2 


2455000.09214 


210.55 


11.67 


2 


2455000.09641 


224.15 


12.23 


2 


2455000.10068 


195.24 


11.97 


2 


2455013.75919 


74.13 


8.47 


2 


2455013.76487 


81.76 


7.53 


2 


2455013.77030 


55.46 


8.11 


2 


2455013.77573 


44.87 


7.60 


2 


2455013.78114 


53.93 


7.71 


2 


2455013.78703 


61.60 


7.50 


2 


2455013.79245 


42.35 


6.83 


2 


2455013.79787 


36.87 


7.31 


2 


2455013.80330 


30.06 


7.10 


2 


2455013.80871 


25.86 


6.95 


2 


2455013.81412 


29.68 


7.44 


2 


2455013.81955 


11.37 


6.86 


2 


2455013.82496 


5.68 


6.69 


2 


2455013.83039 


15.44 


6.63 


2 


2455013.83580 


13.27 


5.85 


2 


2455013.84121 


3.74 


5.69 


2 


2455013.84664 


12.79 


6.91 


2 


2455013.85206 


5.91 


6.67 


2 


2455013.85748 


-4.64 


6.31 


2 


2455013.86290 


0.95 


6.07 


2 


2455013.86832 


-10.16 


6.26 


2 


2455013.87374 


-3.95 


5.32 


2 


2455013.88433 


-7.40 


5.70 


2 


2455013.88976 


-1.91 


6.11 


2 


2455013.89518 


-22.65 


5.91 


2 


2455013.90060 


5.42 


6.32 


2 




^ 1 A 
—J. 14 


^ 9*7 
J.o / 


z 


2455013.91145 


-13.89 


6.19 


2 


2455013.92034 


-14.44 


5.90 


2 


2455013.92575 


-9.34 


5.80 


2 


2455013.93118 


2.95 


5.79 


2 


2455013.93660 


-37.18 


5.53 


2 


2455013.94201 


-22.51 


5.80 


2 


2455013.94744 


-35.31 


5.01 


2 


2455013.95286 


-11.12 


6.17 


2 


2455013.95829 


5.23 


5.75 


2 


2455013.96372 


-59.30 


5.56 


2 


2455013.96913 


-63.07 


6.04 


2 


2455013.97456 


-55.50 


5.88 


2 


2455013.98000 


-64.15 


5.92 


2 


2455013.98542 


-60.13 


5.24 


2 


2455013.99086 


-55.52 


4.96 


2 


2455013.99627 


-69.24 


4.86 


2 



Outlandish orbit of HAT-P-7b 



7 



TAB LEI — Continued 



HJD 


RV [m s 1 ] 


Error [m s 1 ] 


Spec." 


2455014.00486 


-76.79 


5.51 


2 


2455(114 01029 


-80.53 


5.47 


2 




—73.56 


5.96 


2 


2455014 021 14 


—90.86 


5.19 


2 


2455014 02656 


—81.24 


6.44 


2 


2455014 03199 


—109.25 


5.29 


2 


2455014 03742 


—100.31 


5.77 


2 


2455014.04285 


-107.39 


6.12 


2 


2455014.04828 


-107.15 


5.91 


2 


2455014.06427 


-114.12 


5.32 


2 


2455014.06971 


-110.83 


5.39 


2 


2455014.07514 


-114.67 


5.97 


2 


2455014.11226 


-146.78 


5.93 


2 


2455014.11769 


-137.88 


5.84 


2 


2455014.12312 


-143.30 


6.07 


2 


2455014.12855 


-136.48 


5.96 


2 


2455014.13397 


-111.90 


7.02 


2 



NOTE. — The RV was measured relative to an arbitrary template spectrum specific to each spectrograph; only the differences among the RVs from a single spectrograph are 
significant. The uncertainty given in Column 3 is the internal error only and does not account for any possible "stellar jitter." 
a (1) Keck/HIRES, (2) Subaru/HDS. 



TABLE 2 
Model Parameters for HAT-P-7B 



Parameter 


Value 


Orbital period, P [d] 


2.2047304 ± 0.0000024 


Midtransit time [HJD] 


2, 454. 73 1.67929 ±0.00043 


Transit duration (first to fourth contact) [hr] 


4.006 ±0.064 


Transit ingress or egress duration [hr] 


474+0-061 


Planet-to-star radius ratio, R p /Rt 
Orbital inclination, (' [deg] 
Scaled semimajor axis, a/R* 
Transit impact parameter 


0.0834* 
80.8 + f| 

0.618_ n |4 „ 



Velocity semiamplitude, K [m s~'] 


211.8±2.6 


Upper limit on eccentricity (99.73% conf.) 


0.039 


ecosoj 


-0.0019 ±0.0077 


esinw 


0.0037 ±0.0124 


Velocity offset, Keck/HIRES [m s" 1 ] 


-51.2±3.6 


Velocity offset, Subaru/HDS [m s -1 ] 


-4.8 ±2.5 


Constant radial acceleration 7 [m s yr ] 


21.5±2.6 



Projected stellar rotation rate, v sin 4 [km s '] 4 -9^o 9 

Projected spin-orbit angle, A [deg] 182.5 ± 9.4 



